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Abstract. Following the tracks of |Malbet, Yu, fc Shao (1995)| on dark hole algorithms, we 
present analytical methods to measure and correct the speckle noise behind an ideal coron- 
agraph. We show that, in a low aberration regime, wavefront sensing can be accomplished with 
only three images, the next image being fully corrected (no iterative process needed). The only 
hardware required is the coronagraph deformable mirror and an imaging detector in the focal 
plane, thus there are no non-common path errors to correct. Our first method, speckle field 
nulling, is a fast FFT-based algorithm requiring the deformable mirror influence functions to 
have identical shapes. Our second method, speckle energy minimization is more general and 
based on matrix inversion. Numerical simulations show that these methods can improve the 
contrast by several orders of magnitude. 
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systems 



1. Introduction 

By speckle nulling, we refer to any technique whose goal is the reduction of speckle noise 
down to a level compatible with exoplanet detection. In the focal plane, speckle noise 
appears as a background made of a collection of star-like spots, the speckles. Speckles 
doubly hamper planet detection: they worsen the photon noise by increasing the back- 
ground, and they introduce confusion because speckles and planets look alike. 

Speckle noise arises from the scattering of star light by optical defects acting either as 
phase gratings (e.g. polishing errors), or as amplitude gratings (e.g. inhomogeneous reflec- 
tivity). Speckle nulling is achieved by creating destructive interferences at the locations 
of speckles thanks to a deformable mirror (DM for short; see Trauger, Burrows, Gordon, 
et al. 2004). In terms of wave propagation, optical defects cause wavefront aberrations 
that the DM compensate for. 

Speckle nulling belongs to the broader topic of wavefront sensing and control. As 



explained by Malbet, Yu, & Shao (1995) in their pioneering work, it differs from con- 
ventional adaptive optics (AO) in the sense that its goal is not to maximize the Strehl 
ratio (i.e. minimize the standard deviation of phase errors), but to minimize the light 
scattered in a given zone of the image plane, dubbed dark hole, where planet detection 
becomes possible. Moreover, conventional AO measures the wavefront with a separate 
sensing channel, most of the time at a wavelength shorter than the wavelength used 
for science operation. For high-contrast imaging, it is highly desirable to measure the 
wavefront directly with the science image in order to avoid non-common path errors. 
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In their paper, Malbet, Yu, & Shao (1995) discuss only wavefront control assuming 
prior knowledge of the wavefront thanks to some phase retrieval techniques. We provide 
here a full approach to wavefront sensing and control. However, we narrow the discussion 
to the case where aberrations are small, i.e. the problem remains linear. Non-linearity 
could be dealt with by iterating (see Give'on, Kasdin, & Vanderbei 2005) or other meth- 
ods could be used in the first place to get to the wavefront quality needed here. 

This paper is organized as follows: in S|21we expose the wavefront sensing and control 
theory, in jj^jwe present one- and two-dimensional simulations with white and colored 
speckle noise, and we conclude in 21 



2. Speckle nulling theory 

We present here a simplified theory in which diffraction by the pupil is not included. 



For a more complete exposition of the theory, see|Borde fc Traub (2005) 



2.1. Wavefront sensing 

Let E(u, v) — e l ^ w,v ^ « l+icf>(u, v) be the normalized electric field due to the star alone in 
a pupil plane before the coronagraph. Here, is a complex function accounting for small 
departures from a plane wave, i.e. phase and amplitude aberrations. In order to assess the 
performance of speckle nulling alone, we assume an ideal coronagraph which perfectly 
removes the star: in a reimaged pupil after the coronagraph, the field is E'(u,v) — 
i<j)(u,v). On the detector, the image-plane field is then E'(a,(3) = i(j){a,j3), where hats 
denote Fourier transforms. Because the detector measures the field intensity, I(a, (3) = 
\4>(ct, /3)| 2 , complex aberrations cannot be obtained from a single image. Fortunately, in 
an AO system, the DM provides a natural means to acquire the missing information by 
modifying the wavefront phase, which leads to a field intensity variation (this is a form 
of phase diversity). 

It is easy to show that complex aberrations can be measured from a set of three images, 
numbered here 0, 1, and 2. Let ip be the phase change introduced by the DM. Although 
this is not mandatory, we assume for the sake of simplicity that the DM is originally flat 
("00 — 0). The intensity takes the successive values 



/l = i0 + ^|2 ^ 4 )~{* % ~ T ;\~\ 2 (2- 1 ) 
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Provided A = (ipi) ip% — V 7 ! (^2) 7^ 0, we derive for the complex aberrations 
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The condition A ^ reflects that (|2.1I) should not contain twice the same equation, i.e. 
the DM should be driven so that, between images 1 and 2, the intensity is modified in 
a measurable way. In practice, ip\ is computed to modify the intensity in every point by 
an amount comparable to the intensity in that point. Then, 02 is computed to maximize 
A in every point, so as to maximize the measurement signal-to-noise ratio. 

2.2. Wavefront control 

With the knowledge of the complex aberrations, speckles would instantly go away if we 
could make equal to —0. However, 
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(a) the DM cannot take an arbitrary shape as it is made of a facesheet supported 
by a finite number of actuators (see Trauger, Moody, Gordon, et al. 2003 for a detailed 
description), and 

(&) the DM introduces a real, not complex, phase change. 
Point (a) implies that aberrations can only be corrected up to a maximum spatial fre- 
quency set by the density of actuators (with respect to the wavelength), which means 
that the dark hole cannot be arbitrarily large. Although point (6) seems to imply that 
only phase aberrations could be corrected, it really means - as we will show later on - 
that phase and amplitude aberrations can both be corrected, but at the cost of half the 
dark hole area. 

We model the DM shape as the linear superposition of two-dimensional influence func- 
tions. Each of these functions describes the DM surface deformation in response to the 
actuation of a given actuator. The phase change introduced by a iVx TV-element DM is 
then 

JV N 

i>{u,v) = ^2^2a u fki(u,v), (2.3) 

k=0 1=0 

where the a,ki (the actuator strokes) are the unknowns in the wavefront control problem. 

2.2.1. Speckle field nulling 

As a first approach to the ajy calculation, let us examine under what conditions E' 
could be made zero in the simplified case where all influence functions have an identical 
shape, i.e. fki{u, v) = f(u — kj,v — 1 4-), where d is the actuator spacing. We have 

AT-liV-l 2( ff\ 

&a,0+&a I /J)=O <=> ££a H e-^W)=-§^. (2.4) 

k=0 1=0 f Pj 

We recognize in the left-hand size of (|2.4|) a truncated Fourier series. By setting the 
strokes equal to the Fourier coefficients of —<j>/f, we minimize the total energy of E' / f in 
the dark hole. This is close, but not quite the same thing as minimizing the total speckle 
energy in the dark hole, which would be the total energy of E' (more on that in i !2.2.2(l . 
With the Fourier solution to l|2.4|) . we learn that: 

(a) because the Fourier series is truncated (as actuators are in limited number), there 
can be no exact solution to l|2.4|l and no infinitely deep dark hole, unless the truncated 
part of the series happen to be equal to zero, i.e. aberrations are band-limited; 

(b) because Fourier series are periodic, the dark hole will have an extension limited to 
that period. Precisely, it will be a centered square with a maximum size of ±A = ±y-^, 
where D is the pupil's diameter. In other words, a dark hole created with a TV x 7V-DM 
will be at most NxN element of resolution wide, which is consistent with the sampling 
theorem. 

The solution by Fourier series is attractive because it can be coded with FFTs and is 
computationally very fast. However, it has two drawbacks: (1) it cannot handle influence 
function variations across the DM, and (2) it does not exactly minimize the speckle field 
energy in the dark hole which would seem desirable. 

2.2.2. Speckle energy minimization 

As a second approach, let us now compute the au as the result of the minimization of 
the total speckle energy in the dark hole, defined as 

£ = < + ^,0 + ^ > = J J ($+$)($ + $)*. (2.5) 
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Given that dip/daki = fkl, the total energy is minimized when, V(fc, I), 

ni , N-lN-l 

at 



= J] Jh i » 1 !1(</«/h>)=-8(<^>) 1 (2.6) 

flfcZ n=0 m=0 

where 5ft stands for the real part. 

We recognize in ()2.6(l a linear system that can be solved in matrix format: 

FA = $ & A = F' 1 ®, (2.7) 

where F is a TV 2 x iV 2 -matrix, and A and $ are iV 2 -vectors. 

It happens that if influence functions have the property of being separable, that is to 
say if fki{u,v) — gk(u)gi(v), then 1|2.7|) can be rewritten as 

GAG = $^A = G- 1 <Z'G- 1 1 (2.8) 

where G, A and $ are now JVx TV-matrices. This has the advantage of requiring much less 
computing power. Box-shape functions and bidimcnsional Gaussians are two examples of 
fully separable functions. If instrumental influence functions are not fully separable but 
the wavefront is corrected as if they were, a degradation of the dark hole performance is 
to be expected (see A3. 21 for an example). 

When phase and amplitude aberrations are simultaneously present, (f> is complex, and 
the careless use of any of the two algorithms would lead to complex values for the au ■ The 
key to correct amplitude aberrations with real actuator strokes is to treat them as phase 
aberrations. This can be done by deriving from <fi a function with Hermitian symmetry 
and therefore a real Fourier transform. However, that Hermitian function can be made 
equal to <j> only in half the dark hole, therefore half the dark hole area is sacrificed. 



3. Simulations 

3.1. White speckle noise 
To start with, we simulate our speckle nulling algorithms in one-dimension with white 
speckle noise (Fig. QJ. The standard deviation of phase aberrations is set to A/1000 so 



that the linear approximation of i )2.1l would remain valid. Refer to Borde & Traub (2005) 
for simulations with amplitude aberrations as well. The DM features 64 actuators and 
top- hat influence functions. As shown in the bottom panel of Fig. ^ energy minimization 
yields a dark hole slightly deeper than field nulling: 1.4 x 10 -11 vs. 5.8 x 10~ n (inten- 
sities are scaled with respect to the peak of the star Airy function). This is always the 
case by typically a factor of two. It is possible to obtain deeper dark holes with energy 
minimization but at the cost of smaller angular extensions. For instance, if one accepts 
to reduce the dark hole size from 64 to 44 resolution elements, the dark hole floor drops 
down to 2.7 x 10~ 15 . 

3.2. Colored speckle noise 

Figure[3is the result of a two-dimensional simulation including realistic phase aberrations 
from the surface map of an 8.2-m mirror. In this map, the power spectral density of 
surface errors decreases as the third power of the spatial frequency above 0.35 m . 
The DM features 64x64 actuators, and influence functions are as measured on the High- 
Contrast Imaging Testbed (HCIT; Trauger, Burrows, Gordon, et al. 2004) at NASA's Jet 
Propulsion Laboratory. These simulations does not include quantum or read-out noise 
and assume a perfect DM. With field nulling (as shown in Fig. |2J, the average dark hole 
floor is 5.9 x 10~ 12 . Energy minimization assuming separable influence functions yields 
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Figure 1. One-dimensional speckle nulling simulation. Top: Phase-aberrated wavefront (dotted 
line) and negative of DM shape (thick solid line). Bottom: Uncorrected coronagraphic image 
(dotted line), and corrected coronagraphic image with speckle field nulling (thin solid line) and 
with speckle energy minimization (thick solid line). 

a higher average floor of 7.1 x 10 -11 , because of a 5 % mismatch between the real HCIT 
functions and their separable approximations. 



4. Conclusion 

In this paper, we presented a full theory for wavefront sensing and control using sci- 
ence camera images (no separate wavefront sensing channel), assuming that phase and 
amplitude aberrations remain small. Sensing is performed with only three images, and 
relies on the deformable mirror to introduce the needed phase diversity. 

Since our theory is general, it can adapted to most coronagraph designs. A good exam- 



ple is the application to shaped pupil coronagraphs by Give'on, Kasdin, & Vanderbei (2005) 
We are currently working on the adaptation to coronagraphs with image-plane band- 
limited masks (Kuchner & Traub 2002) to start experimenting with the HCIT. In paral- 
lel, we plan to develop a polychromatic theory to accommodate the spectral bandwidth 
necessary for realistic planet detection and spectroscopy. 
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Uncorrected coronagraphic image Dark hole (64 x 64 actuators) 




Figure 2. Two-dimensional speckle nulling simulation. Bottom left: phase map of the first 8.2-m 
primary mirrors of ESO's Very Large Telescope. Top left: uncorrected coronagraphic image 
obtained with this phase map scaled to a standard deviation of A/1000. Top right: corrected 
image obtained with the field nulling algorithm ( H2.'2.1t . Bottom right: average cut through the 
dark hole. 
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Discussion 

Trauger: Your model assumes that coronagraph optics downstream from the DM do not 
introduce additional phase and amplitude errors. However, these errors will necessarily 
exist and affect the image. Could you comment on that? 

Borde: We have indeed assumed an optically perfect coronagraph. I guess that the 
solution to this problem will be to iterate the algorithm, so that unknown errors will be 
taken care of in the process. In addition, I should say that the effect of the coronagraph 
is not only to remove the star PSF as in our simplistic model, but is also to modify the 
DM influence functions in a way that depends on their positions in the pupil. This can 
easily be (and will be) incorporated in the energy minimization algorithm, but cannot 
be dealt with in the field nulling algorithm. 



